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SYSTEM AND METHOD FOR DEFINING T-SPLINE AND 
T-NURCC SURFACES USING LOCAL REFINEMENTS 



The patent application claims priority to U.S. Provisional Patent Application No. 
5 60/458,23 1 , entitled "T-splines" and filed on March 26, 2003. 

FIELD OF THE INVENTION 

The present invention relates generally to defining a modeled surface that is capable 
of expressing local refinements. 
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BACKGROUND 

Surface modeling is a fundamental task in computer graphics, computer aided design 
(CAD), computer aided geometric design (CAGD), and computer animation. Surfaces are 
typically approximated by a mesh of polygons or piecewise polynomial patches. Because of 

1 5 their simplicity, polygonal meshes are a popular way to approximate surfaces and are 

adequate for many applications. A drawback of polygonal meshes is that they are inherently 
faceted, and a large number of polygons may be needed to make the facets small enough to 
satisfy the demands of the application. 

In contrast, a popular smooth surface representation is tensor product B-spline 

20 surfaces. The control points for a B-spline surface are required to be topologically arranged 
in a rectangular grid. B-spline surfaces are comprised of several parametric surface patches 
that can be represented in Bezier form. A significant advantage of B-spline surfaces is that 
each of the constituent surface patches are automatically C n_1 with its neighboring patches, 
where n is the degree of the basis functions. 

25 The parameter values at which the constituent surface patches begin and end are 

called knots, and a non-decreasing sequence ofknots is called a knot vector. A B-spline 
surface definition includes two knot vectors, one for each parameter of the parametric 
equation. In a uniform B-spline surface, the difference between each pair of knots in a given 
knot vector is constant. A non-uniform B-spline does not have such a constraint. A rational 

30 B-spline surface is one for which the control points are assigned weights. Weights provide 
additional shape control, allow the introduction of sharp creases in the model, and empower 
the B-spline surface to represent quadric surfaces. A non-uniform rational B-spline surface is 
referred to by the acronym NURBS surface. 



2 

Refinement of a NURBS control grid is accomplished through a procedure called knot 
insertion whereby one or more new knot values are inserted into a knot vector. Knot 
insertion does not alter the surface, but it provides more control points with which the 
designer can manipulate the model. Because NURBS control points must lie topologically in 
5 a rectangular grid, knot insertion causes one or more entire rows of control points to be added 
to the control grid. Cartesian coordinates for those new rows of control points along with 
some of the neighboring control points can be found such that the refinement operation does 
not alter the surface. Unfortunately, adding an entire row of control points for each desired 
new control point simply to satisfy the rectangular grid topology of NURBS surfaces 

10 increases the modeling complexity of the surface dramatically. 

Surface refinement is a valuable operation that has several uses. First of all, 
refinement allows a designer to insert additional control points into a surface region where 
more detail is called for in the model. For example, more control points would be needed 
using NURBS to model a human ear than to model a cheek. Second, each time a control grid 

1 5 is refined, the control grid itself becomes an increasingly better approximation of the NURBS 
surface that it defines. Thus, by performing repeated refinements, the control grid can 
become an arbitrarily accurate representation of the surface, which is suitable for rendering. 
Third, knot insertion can be used to compute the control points of the Bezier patches that 
comprise a NURBS. Fourth, sharp features can be added to NURBS surfaces using knot 

20 insertion. NURBS surfaces are generally C n_1 continuous, where n is the degree of the 

surface. However, if r identical knots exist in a knot vector, the surface is C n " r at that knot 
value. Thus, a triple knot in a cubic NURBS surface permits a C° crease. Uniform B-spline 
surfaces do not allow arbitrary knot insertion, for as soon as a single knot is inserted, the 
surface becomes a non-uniform B-spline surface. The only knot insertion possible for a 

25 uniform B-spline surface (if it is to remain uniform after the knot insertion) is to 

simultaneously insert a knot midway between each pair of neighboring knots, a process called 
knot doubling. Since uniform B-splines cannot support two or more identical knots, uniform 
B-spline surfaces are always C n_1 continuous and it is not possible to impose a C° crease. 

Several algorithms exist for knot insertion. The Oslo algorithm computes so-called 

30 discrete B-splines to define the B-spline basis transformation from a refined space of splines 
to a subspace. Another method of knot insertion is Boehm's algorithm which works directly 
with the B-spline coefficients. In addition, mathematical insights like the blossoming 
principle for knot insertion have been developed. 
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Knot removal, the inverse operation of knot insertion, can also be performed on 
NURBS. By comparison, knot insertion does not modify the surface, while knot removal in 
general does alter the surface. Therefore, knot removal usually involves approximation. One 
motivation for knot removal is data reduction. The problem is to minimize the number of 
5 knots subject to an error tolerance. Another application of knot removal is shape fairing. The 
continuity order can be increased by removing knots. 

The requirement that all control points for a NURBS surface must topologically be 
positioned in a rectangular grid is a serious shortcoming. There are at least three reasons why 
this is a problem. First, surfaces of arbitrary topology can be represented using NURBS 

10 surfaces only by partitioning the model into a collection of individual NURBS patches. 

Adjacent patches are then explicitly stitched together using geometric continuity conditions. 
FIG. 2 shows a hand model comprised of seven NURBS surfaces. The small rectangular area 
is magnified in FIG. 3 to show a hole where neighboring NURBS surfaces do not match 
exactly. The presence of such gaps places a burden on modelers who potentially must repair 

1 5 a widened gap whenever the model is deformed. 

A second reason why this constraint is so serious is because it typically means that a 
large number of NURBS control points serve no purpose other than to satisfy topological 
constraints. They carry no significant geometric information. Superfluous control points are 
a serious nuisance for designers, not merely because they require the designer to deal with 

20 more data, but also because they can introduce unwanted ripples in the surface. FIG. 1 shows 
a NURBS head model. Designers can waste dozens of hours on models, such as FIG. 1, in 
tweaking the NURBS control points while attempting to remove unwanted ripples. The 
darker NURBS control points are superfluous (as will be later shown in the present invention) 
arid have been added as the result of including control points desired by the modeler that 

25 resulted in additional rows of control points. 

A third reason why the rectangular grid constraint is a problem is that it is not possible 
to insert a single control point into a NURBS control grid, but rather, an entire row of control 
points must be added each time a single knot insertion is performed. This means that true 
local refinement of NURBS surfaces is not possible, since the insertion of a single control 

30 point demands the insertion of an entire row into the control grid. 

In order to further illustrate the complexities of solving the problems associated with 
NURBS, one attempt at solving this problem will now be addressed. In August 2001, Almaz 
Bakenov completed a master's thesis at Brigham Young University entitled, "T-splines: 
Tensor Product B-spline Surfaces with T- Junctions". Thomas W. Sederberg, the inventor of 
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the present invention described later in the detailed description, was the advisor for Almaz 
Bakenov's thesis. The Bakenov thesis represents an early attempt at resolving the problems 
associated with NURBS. Indeed, Thomas W. Sederberg wrestled with this problem from 
1999 to 2003. Thus, we have the slightly confusing situation in which there are two 
5 significantly different concepts called T-splines: the concept later described in the present 
invention, and the notion of T-splines described in the Bakenov thesis. 

Particularly, the notion of T-splines presented in the Bakenov thesis contains 
limitations that render it almost useless. Most significantly, the T-splines described in the 
Bakenov thesis are not capable of solving the problems associated with NURBS without the 

1 0 solution of a potentially huge system of linear equations, for which a solution does not exist 
in many cases. For these reasons, the notion of T-splines contained in the Bakenov thesis is 
of little practical value. 

Several surface formulations have been proposed that do not suffer from NURBS' 
topological constraint of requiring all control points to lie in a rectangular grid. The 

15 arbitrary-topology surface formulation that is most pertinent to the present discussion is 

called a subdivision surface. Catmull-Clark subdivision surfaces are, in fact, generalizations 
of B-spline surfaces. When the control points of a Catmull-Clark surface happen to lie 
topologically in a rectangular grid, the surface degenerates to a bicubic uniform B-spline 
surface. In like manner, Doo-Sabin subdivision surfaces generalize biquadratic uniform B- 

20 spline surfaces to control grids of arbitrary topology. 

Refinement of a Catmull-Clark surface is based on the notion of knot doubling for 
uniform bicubic B-spline surfaces, with special rules introduced to handle non-four-sided 
faces and non- valence-four vertices. With each refinement step, the number of faces grows 
by a factor of four. Repeated refinement causes the control grid to approximate the limit 

25 surface ever more closely. 

A problem with Catmull-Clark surfaces is that, since they are based on uniform B- 
spline surfaces, it is not possible to represent sharp C° features. However, some methods 
have been proposed for defining sharp features, edges, creases, and corners into subdivision 
surfaces by altering the refinement rules in the neighborhood of such features. 

30 Another more serious problem is that local refinement of Catmull-Clark surfaces is 

not possible. In fact, where the simplest refinement operation for a NURBS surface is to 
insert a single row of control points, only a global refinement operation is defined for 
Catmull-Clark surfaces, and the number of control points grows by roughly a factor of four 
with each refinement step. 
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SUMMARY OF THE INVENTION 

A system and method is provided for creating a bi-cubic spline surface in a computing 
5 environment. One operation in the method is defining a control mesh with a substantially 
rectangular structure. A further operation is inferring from the control mesh the tensor 
product B-spline basis functions for each control point. The surface can then be computed 
based on the basis functions and the control mesh. 

Additional features and advantages of the invention will be apparent from the detailed 
10 description which follows, taken in conjunction with the accompanying drawings, which 
together illustrate, by way of example, features of the invention. 



BRIEF DESCRIPTION OF THE DRAWINGS 

1 5 FIG. 1 is a perspective view of a NURBS model of a human head; 

FIG. 2 illustrates a NURBS model of a human hand; 

FIG. 3 depicts a hole where neighboring surfaces do not exactly match in the NURBS 
model of the hand in FIG. 2; 

FIG. 4 depicts the two surfaces in the NURBS model human hand of FIG. 2 that have 
20 been corrected using a T-spline using an embodiment of the present invention; 

FIG. 5 illustrates a head modeled as a NURBS surface and a head modeled as a T- 
spline surface in an embodiment of the invention; 

FIG. 6 illustrates an embodiment of a Catmull-Clark mesh that is refined using T- 
NURCC local refinements; 
25 FIG. 7 is a flow chart illustrating an embodiment of a method for defining a bi-cubic 

spline surface in a computing environment; 

FIG. 8 depicts an example of a cubic B-spline curve with a knot vector; 

FIGS. 9a and 9b illustrate a region of a NURBS control mesh labeled with knot 
intervals; 

30 FIG. 10 shows an embodiment of a pre-image of a portion of a T-mesh in (s,t) 

parameter space; 

FIG. 1 1 illustrates knot lines for a basis function Bj(s,t); 

FIG. 12 illustrates a sample refinement of Bi(s,t) in an embodiment of the present 
invention; 
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FIG. 13 illustrates an embodiment of a nested sequence of T-spline spaces; 
FIGS. 14a - 14f illustrate a local refinement example in an embodiment of the 
invention; 

FIG. 1 5 illustrates the insertion of a dart in a T-mesh in an embodiment of the 
5 invention; 

FIG. 16 illustrates semi-standard T-splines in an embodiment of the invention; 

FIG. 17 is a flowchart illustrating an embodiment of a method defining bicubic spline 
surfaces that provides local refinement to control meshes in a computing environment; 

FIG. 1 8 illustrates Bezier domains in a pre-image of a T-mesh in an embodiment of 
10 the present invention; 

FIGS. 19a and 19b illustrate the creation of Bezier control points using local knot 
insertion; 

FIG. 20 illustrates the merging of two B-splines; 

FIG. 21 illustrates merging two B-splines using cubic NURSSes and the problems 
15 created by such a merger; 

FIG. 22 depicts the merging of two B-splines using T-splines in an embodiment of the 
present invention; 

FIG. 23 depicts the results of the merging of two surfaces using B-splines ? C° T- 
splines and C l T-splines in an embodiment of the present invention; 
20 FIG. 24 illustrates local refinement using T-splines about a valence-four control point; 

FIG. 25 illustrates local refinement using T-splines about a n-valence control point; 

and 

FIG. 26 illustrates a T-NURCC showing the influence of parameter p where the limit 
surface of the control grid on the top left isp = 0.1 and the top right is p = 0.5 and the bottom 
25 right is p = 0.9. 

DETAILED DESCRIPTION 

Reference will now be made to the exemplary embodiments illustrated in the 
30 drawings, and specific language will be used herein to describe the same. It will nevertheless 
be understood that no limitation of the scope of the invention is thereby intended. Alterations 
and further modifications of the inventive features illustrated herein, and additional 
applications of the principles of the inventions as illustrated herein, which would occur to one 
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skilled in the relevant art and having possession of this disclosure, are to be considered within 
the scope of the invention. 

The present system and method of the invention provide a generalization and 
improvement of non-uniform B-spline surfaces (NURBS) called T-splines. A major problem 
5 with NURBS surfaces is that, because NURBS control points must lie topologically in a 
rectangular grid, it is often the case that a large percentage of the NURBS control points 
serve no purpose other than to satisfy the rectangular grid topology. We refer to such control 
points as superfluous, because they contain no significant geometric information. T-splines 
are a generalization of NURBS surfaces that are capable of minimizing the number of 

10 superfluous control points. In addition, T-splines can insert a control point into the control 
mesh such that a bi-cubic spline surface is not geometrically altered. 

One difference between a T-mesh (or a T-spline control mesh) and a NURBS control 
mesh is that T-splines allow a row of control points to terminate.. Accordingly, a row can 
include many points or just a single point. The final control point in a partial row is called a 

1 5 T-junction. The T-junctions can be seen in FIG. 13b and c as compared to the rectangular 
grid in FIG. 13a and d. FIG. 5b shows another example of a T-spline in which the 
superfluous control points in a NURBS head model of FIG. 5a are removed or avoided when 
using a T-spline. A T-spline surface model is geometrically equivalent to a NURBS model, 
yet a T-spline representation can often be made using substantially fewer control points than 

20 an equivalent NURBS model. 

T-spline control grids permit T-junctions, so lines of control points need not traverse 
the entire control grid as with a NURBS. In other words, T-junctions allow T-splines to be 
locally refineable and control points can be inserted into the control grid without propagating 
an entire row or column of control points. 

25 T-splines support many valuable operations within a consistent framework, such as 

local refinement, and the merging of several B-spline surfaces that have different knot vectors 
into a single gap-free model. This detailed description focuses on T-splines of degree three, 
which are C 2 in the absence of multiple knots. However, the use of T-splines extends to any 
degree. T-NURCCs (Non-Uniform Rational Catmull-Clark Surfaces with T-junctions) are a 

30 superset of both T-splines and Catmull-Clark surfaces. Thus, T-NURCCs will also be 

discussed in this disclosure. In addition, a modeling program for T-NURCCs can handle any 
NURBS or Catmull-Clark model as special cases. 

T-NURCCs enable true local refinement of a Catmull-Clark-type control grid. This 
means individual control points can be inserted only where they are needed to provide 
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additional control, or to create a smoother tessellation, and such insertions do not alter the 
geometry of the surface or the limit surface. T-NURCCs use stationary refinement rules and 
are C 2 except at extraordinary points and features. The detailed description also presents a 
locally refineable subdivision surface called T-NURCCs (Non-Uniform Rational Catmull- 
5 Clark surfaces with T-junctions). In T-NURCCs, faces adjacent to an extraordinary point can 
be refined without propagating the refinement, and faces in highly curved regions can also be 
refined locally. As in T-splines, individual control points can also be inserted into a T- 
NURCC to provide finer control over details. T-NURCCs are a generalization of both 
Catmull-Clark surfaces and NURBS. 

10 FIG. 6 shows how a T-NURCC local refinement enables a T-NURCC tessellation to 

be far more economical than a globally-refined Catmull-Clark surface. The T-NURCC 
version of the tetrahedral shape in FIG. 6 has 2496 faces. A globally refined Catmull-Clark 
surface for the tetrahedral shape needs 393,216 faces to achieve the same precision. 

Returning again to FIG. 2, the figure shows a hand model comprised of seven B- 

1 5 spline surfaces. The small rectangular area is blown up in FIG. 3 ,to magnify a hole where 
neighboring B-spline surfaces do not match exactly. The presence of such gaps places a 
burden on modelers, who may need to repair a widened gap whenever the model is deformed. 
FIG. 3 shows the gap created using NURBS and FIG. 4 shows the model after being 
converted into a gap-free T-spline, thereby eliminating the need for repair. T-splines and T- 

20 NURCCs can thus imbue models comprised of several NURBS surfaces with the same air- 
tightness that Catmull-Clark surfaces extend to uniform cubic B-spline based models. T- 
splines are an enhancement of NURBS surfaces that allow the presence of T-junction control 
points. 

FIG. 7 illustrates the fundamental nature of the present invention and provides an 
25 overview thereof. A method is provided for defining a bi-cubic spline surface in a computing 
environment. One operation in the method is creating a control mesh with a substantially 
rectangular structure as in block 150. Another operation is inferring from the control mesh 
the tensor product B-spline basis functions for each control point as in block 152. A further 
operation is computing the surface based on the basis functions and the control mesh as in 
30 block 154. The basis function can also be determined for each control point using one non- 
hierarchical set of rules. 

Refinement of a control mesh is the process of adding control points to the control 
mesh without altering the surface defined by the control mesh. As discussed, the simplest 
form of refinement for NURBS surfaces requires the addition of an entire row of control 
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points to the control grid. In contrast, the present invention enables a modeler to refine a T- 
spline control grid with a single control point. In simple terms, the surface as it existed before 
the addition of the control point is unaltered. Prior known surfaces such as NURBS do not 
support the addition of a single new control point. This problem will be illustrated in further 
5 detail later. The present invention has been discussed in general terms up to this point but a 
more technical description of the invention will now follow. 

Knot Intervals 

1 0 A knot interval is a non-negative number assigned to each edge of a T-spline control 

grid for the purpose of conveying knot information. In the cubic B-spline curve shown in 
FIG. 8, the dj values that label each edge of the control polygon are knot intervals. Each knot 
interval is the difference between two consecutive knots in the knot vector. For a non- 
periodic curve, end-condition knot intervals are assigned to "phantom" edges adjacent to each 

1 5 end of the control polygon (in this case d.i and d 5 ). For all but the first and last edges of the 
control polygon, the knot interval of each edge is the parameter length of the curve segment 
to which the edge maps. Any constant can be added to the knots in a knot vector without 
changing the knot intervals for the curve. Thus, if the knot intervals are given and a knot 
vector is desired to be inferred, then a knot origin can be chosen. 

20 Edges of T-spline and T-NURCC control grids are likewise labeled with knot 

intervals. Since T-NURCC control meshes are not rectangular grids, knot intervals enable 
local knot coordinate systems to be imposed on the surface. FIG. 9 shows a regular subgrid of 
a NURCC control grid. We can impose a local knot coordinate system on this region, and 
therewith determine basis functions for the control points, as follows. First, (arbitrarily) 

25 assign P 0 o the local knot coordinates of (J 0 > ^o). The local knot vectors for this regular subgrid 
are then {d} = {0, d 0 , d x , . . . } and {e} = {0, e 0 , e x , . . . } where 

_ / i 

30 and =e_ } = 0. The polar label for control point P,y, with respect to this local knot 

coordinate system, is thus /(5 M ,5 p 3 i+! ;c y _ 1 ,e y ,e /+1 )and the surface defined by this regular 
subgrid is 
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p(j.o=EEv?(*)n}(o 

' J 

where N? (s) are the cubic B-spline basis functions over {d} and the (t) are over {e} . The 
superscript 3 denotes degree, not order. 

5 

T-Splines 

A control grid for a T-spline surface is called a T-mesh. If a T-mesh forms a 
rectangular grid, the T-spline degenerates to a B-spline surface. A T-mesh is basically a 
10 rectangular grid that allows T-junctions. The pre-image of each edge in a T-mesh is a line 
segment of constant s (which we will call an s-edge) or of constant / (which we will call a t- 
edge). A T-junction is a vertex shared by one s-edge and two f-e4ges, or by one f-edge and 
two sledges. 

A knot interval is assigned to each edge in the T-mesh. FIG. 10 shows the pre-image 
15 of a portion of a T-mesh in (s,t) parameter space, where the d* and ej denote the knot intervals. 
Knot intervals are constrained by the relationship that the sum of all knot intervals along one 
side of any face must equal the sum of the knot intervals on the opposing side. For example, 
in FIG. 10 on face Fi, e3+e4 = e6+e7, and on face F 2 , d6+d7 = d9. 

The inventor of the present invention has realized that it is possible to infer a local 
20 knot coordinate system from the knot intervals on a T-mesh. To impose a knot coordinate 
system, a control point can be first chosen whose pre-image will serve as the origin for the 
parameter domain (s,t) = (0,0). For the example in FIG. 10, (s 0 ,to) can be designated to be the 
knot origin. 

Once a knot origin is chosen, an s knot value can be assigned to each vertical edge in 
25 the T-mesh topology, and a t knot value to each horizontal edge in the T-mesh topology. In 
FIG. 10, those knot values are labeled Si and tj. Based on the choice of knot origin, we have 
s 0 = to = 0, si = di, s 2 = di+d 2 , s 3 = dj+d 2 +d 3 , ti = ei, t 2 = ei+e 2 , and so forth. Likewise, each 
control point has knot coordinates. For example, the knot coordinates for Pi are (s 2 , t 2 +e6), 
for P 2 are (ss,t 2 ), and for P3 are (ss,t 2 +e 6 ). 
30 One additional rule for T-meshes is that if a T-junction on one edge of a face can 

legally be connected to a T-junction on an opposing edge of the face (thereby splitting the 
face into two faces), that edge must be included in the T-mesh. Legal means that the sum of 
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knot vectors on opposing sides of each face must always be equal. Thus, a horizontal line 
would need to split face Fi if and only if e3 = e6 and therefore also e 4 =e7. 

The knot coordinate system is used in writing an explicit formula for a T-spline 
surface: 

(1) P(*,0 = (x(s 9 0,y(s,0,z(s,OMs,t)) = £P,*,(J,0 

where Pi = (xj,yi,Zi,Wi) are control points in P 4 whose weights are Wj, and whose Cartesian 
coordinates are (xj/wj,yi/wj,Zj/wj). Likewise, the Cartesian coordinates of points on the 
surface are given by 



£(x i ,y i ,z i )5 / (j , ,0 
10 (2) 



The basis functions Bj(s,t) are given by 

Bi(s,t) = N[Si0 5 Si 1 ,Si2,Si3,Si4](s)N[t i o,tii,ti2,ti3,ti 4 ](t) 

where N[Sio,Sii 5 Si2,Si3,Si 4 ](s) is the cubic B-spline basis function associated with the knot vector 
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(3) Sj = [SiO,Sii,Si2,Si3,Si 4 ] 

and N[tio,tii,ti2,ti3,ti 4 ](t) is associated with the knot vector 

20 (4) ti=[tio,tii,ti2,ti 35 t i4 ] 

as illustrated in FIG. 11. 

The designer is free to adjust the weights wj to obtain additional shape control, as in 
rational B-splines. As explained later, weights also play an important role in our new local 
25 refinement algorithm: A T-spline whose weights are all one before refinement might end up 

n 

with some weights not equal to one after refinement, yet ]T w .B£s 9 1) = 1 . 

i=i 

The T-spline equation is similar to the equation for a tensor-product rational B-spline 
surface. One difference between the T-spline equation and a B-spline equation is in how the 
knot vectors Sj and t x are determined for each basis function Bi(s,t). Knot vectors Si and ti are 
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inferred from the T-mesh neighborhood of Pi. Since reference will be made to the rule 
whereby the knot vectors are inferred, it is formally stated as: 

Rule 1. Knot vectors Si and tj for the basis function of Pi are determined as follows. 
5 (Si2,ti2) are the knot coordinates of Pj. Consider a ray in parameter space R(a) = (Sj2+a,tj2). 
Then Si3 and Sj4 are the s coordinates of the first two s-edges intersected by the ray (not 
including the initial (Sj2,ti2). By s-edge, we mean a vertical line segment of constant s in the 
pre-image of a T-mesh. The other knots in s t and ti are found in like manner. 

Rule 1 can be illustrated by a few examples. The knot vectors for Pi in FIG. 10 are si 
10 = [so,Si,S2 5 S3,s 4 ] and ti= [ti,t2,t2+e 6 ,t4 5 t 5 ]. For P 2 , s 2 = [s 3 ,S4,S5,S6,s 7 ] and t 2 =[to 5 ti,t2,t2+e 6 ,t4]. For 
P3, s 3 = [s 3 ,S4,S5,S7 5 s 8 ] and t 3 = [ti,t2,t2+e 65 t4 5 t 5 ]. Once these knot vectors are determined for each 
basis function, the T-spline is defined using the above-mentioned T-spline equations. 

T-spline Local Refinement 

1 5 An overview of one embodiment of operations of the invention for inserting control 

points to for a T-spline will now be described. It should also be realized that other 
embodiments of this invention exist and that the invention can also be implemented in 
different forms that are computationally equal. The surface and its associated control mesh 
are typically calculated in Cartesian coordinates. However other coordinate systems can be 

20 used such as Polar coordinates, curvilinear coordinates, or one of the many other known 

coordinate systems. Basis function refinement plays an important role in this algorithm and 
is reviewed first. Next, the notion of T-spline spaces is introduced. This concept is used in 
the local refinement algorithm. 

25 

Basis Function Refinement 

If s = [so,Si,s 2 ,S3 5 s 4 ] is a knot vector and s is a knot vector with m knots with s a 
subsequence of S , then N[so 5 Si,s 2 ,S3,s 4 ](s) can be written as a linear combination of the m-4 
B-spline basis functions defined over the substrings of length 5 in S. 
30 We now present all basis function refinement equations for the knot case m=6. These are the 
cases for which a single knot is inserted into s to create S . Equations for m>6 can be found by 
repeated application of these equations. If s= [so 5 k,si,s 2 ,S3,s 4 ] then 

(5) N(s) = (k - s 0 )/(s3 - so) N[so,k,si,S2,s 3 ](s) + N[k, Si,s 2 ,s 3 ,s 4 ](s) 
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If S= [So,Si,k,S 2 ,S 3 ,S 4 ], 

(6) N(s) = (k-s 0 )/(s 3 -s 0 ) N[so,Si,k,s 2 ,s 3 ](s) + (s 4 -k)/(s4-si)N[si J k,S2,S3 > s 4 ](s) 

5 If S= [So,Si,S 2 ,k,S 3 ,S 4 ], 

(7) N(s) = (k-s 0 )/(s 3 -s 0 ) N[s 0 ,si,s 2 ,k,s 3 ](s) +(s 4 -k)/(s 4 -si) N[si,s 2 ,k,s 3 ,s 4 ](s) 

If S = [So>Si 5 S2,S 3 ,k,S 4 ], 

10 

(8) N(s) = N[s 0 ,si,s 2 ,s 3 ,k](s) + (s 4 -k)/(s 4 -Si)N[si 5 s 2 ,s 35 k ? s 4 ](s) 

If k < s Q or k > s 4 , N(s) does not change. 

A T-spline basis function B(s 5 t) can undergo knot insertion in either s or t, thereby 
15 splitting it into two scaled basis functions that sum to the initial function. Further insertion 
into these resultant scaled basis functions yields a set of scaled basis functions that sum to the 
original function. 

For example, FIG. 12a shows the knot vectors for a T-spline basis function Bi, and 
FIG. 12b shows a refinement of the knot vectors in FIG. 12.a. By appropriate application of 
20 (5) — (8), we can obtain 

(9) B x (s,t) = c\B x (s,t) + c?B 2 (s,t) + clB^sJ) + c A x B, (s,t) 

25 T-spline Spaces 

A T-spline space can be defined as the set of all T-splines that have the same T-mesh 
topology, knot intervals, and knot coordinate system. Thus, a T-spline space can be 
represented by the diagram of a pre-image of a T-mesh such as in FIG. 10. Since all T-splines 
in a given T-spline space have the same pre-image, it is proper to speak of the pre-image of a 
30 T-spline space. A T-spline space Si can be described as a subspace of S 2 (denoted S 2 => Si) if 
local refinement of a T-spline in Si will produce a T-spline in S 2 (discussed later in the 
section on local refinement). If Ti is a T-spline, then Ti eSi means that Ti has a control grid 
whose topology and knot intervals are specified by Si. 
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FIG. 13 illustrates a nested sequence of T-spline spaces, that is, S { a S 2 <z S 3 a ... <= 5*,,. 
Given a T-spline P(s,t) e S l9 denote by P the column vector of control points for P(s,t). 
Given a second T-spline P(s y t) e S 2 , such that P(s,t) = P(s,t), denote by P the column 
vector of control points for P(s,t). There exists a unique linear transformation that maps P 
5 into P . We can denote the linear transformation 

(10) M l2 P = P. 

The matrix Mi, 2 is found as follows. P(s,t) is given by (1), and 

10 

(11) P(5,0 = ZPA(*0 

As described in (9), each B_i(s,t) can be written as a linear combination of the Bj(s,t): 

n 

15 (12) BjfcO^Z*/ 5 /*') 

We require that the refined surface is equivalent to the initial surface: P(s,t) = P(s,t). This 

n 

requirement is satisfied if P y = Y, c i P i • Thus ' the element at row j ^ column i of Mi j2 in 

/=i 

(10) is c, j . In this manner, it is possible to find transformation matrices My that maps any T- 

20 spline in Si to an equivalent T-spline in Sj, assuming S i a 5 y . 

The definition of a T-spline subspace S i c Sj means more than simply that the 
preimage of Sj has all of the control points that the preimage of Sj has, because sometimes it 
is not possible to add a given control point to an existing T-mesh- without adding other control 
points as well. The section on local refinement presents insight into why that is, and presents 

25 our local refinement algorithm for T-splines. This, of course, will allow us to compute valid 
superspaces of a given T-spline space. 
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Local Refinement Method 

T-spline local refinement means to insert one or more control points into a T-mesh 
without changing the shape of the T-spline surface. This procedure can also be called local 
knot insertion, since the addition of control points to a T-mesh is accompanied by 
knots inserted into neighboring basis functions. 

This embodiment of the refinement algorithm has two phases: the topology phase, 
and the geometry phase. The topology phase identifies which (if any) control points must be 
inserted in addition to the ones requested. Once all the required new control points are 
identified, the Cartesian coordinates and weights for the refined T-mesh are computed using 
the linear transformation presented in the section on T-spline spaces. We now explain the 
topology phase of the algorithm. 

An important concept in understanding this discussion is to keep in mind how in a T- 
spline, the basis functions and T-mesh are tightly coupled: For every control point there is a 
corresponding a basis function and each basis function's knot vectors are defined by Rule 1 . 
In this invention, the basis functions are temporarily decoupled from the T-mesh. This means 
that during the computation method, the existence of basis functions is permitted that violate 
Rule 1, and control points may temporarily exist to which no basis functions are attached. 

Our discussion distinguishes three possible violations that can occur during the course 
of the refinement algorithm: 

Violation 1. A basis function is missing a knot dictated by Rule 1 for the current T-mesh. 
Violation 2. A basis function has a knot that is not dictated by Rule 1 for the current T-mesh. 
Violation 3. A control point has no basis function associated with it. 

If no violations exist, the T-spline is valid. If violations do exist, the algorithm resolves them 
one by one until no further violations exist. Then a valid superspace has been found. 
The topology phase of our local refinement algorithm consists of these steps: 

1 . Insert all desired control points into the T-mesh. 

2. If any basis function is guilty of Violation 1, perform the necessary knot insertions into 
that basis function. 
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3. If any basis function is guilty of Violation 2, add an appropriate control point into the T- 
mesh. 

4. Repeat Steps 2 and 3 until there are no more violations. Resolving all cases of Violation 1 
and 2 will automatically resolve all cases of Violation 3. 

5 

We illustrate the method with an example. FIG. 14a shows an initial T-mesh into 
which we wish to insert one control point, P 2 . Because the T-mesh in FIG. 14a is valid, there 
are no violations. But if we simply insert P2 into the T-mesh (FIG. 14b) without changing any 
of the basis functions, we introduce several violations. Since P2 has knot coordinates (S3,t2>, 
10 four basis functions become guilty of Violation 1 . Specifically, those centered at (sj,t 2 ), 

(S2>t2), (S4,t2) ? and (s 5 ,t2). To resolve these violations, we must insert a knot at S3 into each of 
those basis functions, as discussed in the section on basis function refinement. The basis 
function centered at (s 2 ,t 2 ) is N[s 0) Si,S2 5 S4,S5](s)N[to,ti,t2 5 t 3 ,t4](t). Inserting a knot s=s 3 into the 
s knot vector of this basis function splits it into two scaled basis functions: 

15 (S3-So)/(S4-So)N[So ? S 1 ,S2,S3 5 S4](s)N[to J t 1 ,t 2 ,t3,t 4 ](t) (FIG. 14c) 

and (s4-S3)/(s4-Si)N[si,S2,S3,S4,s 5 ](s) N[to,ti,t 25 t 3 ,t4](t) (FIG. 14d) as given in (7). 

The basis function (s3-so)/(s4-So)N[so,si f S2,S3,S4](s)N[tb,ti,t2,t 3 ,t4](t) in FIG. 14c satisfies Rule 

1. Likewise, the refinements of the blending functions centered at (si,t 2 ), (s 4 ,t 2 ), 

and (s 5 ,t 2 ) all satisfy Rule 1. However, the t knot vector of basis function d 2 
20 N[si,S2,s 3 >S4 5 s 5 ](s) N[to 3 ti,t2 5 t 3 ,t4](t) shown in FIG. 14d is guilty of Violation 2 because the 

basis function's t knot vector is [to,ti,t 2 ,t3,t4], but Rule 1 does not call for a knot at t 3 . This 

problem cannot be remedied by refining this basis function, but an additional control point 

can be added into the T-mesh to resolve the violation. 

The needed control point is P 3 in FIG. 14e. Inserting that control point fixes the case 
25 of Violation 2, but it creates a new case of Violation 1. As shown in FIG. 14f, the blending 

function centered at (s 2 ,t 3 ) has an s knot vector that does not include s 3 as required by Rule 1 . 

Inserting s 3 into that knot vector fixes the problem, and there are no further violations of Rule 

1. 

This method will terminate, because the basis function refinements and control point 
30 insertions must involve knot values that initially exist in the T-mesh, or that were added in 
Step 1 . In the worst case, the algorithm would extend all partial rows of control points to 
cross the entire surface. In practice, the algorithm typically requires few if any additional new 
control points beyond the ones the user wants to insert. 
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Local knot insertion is useful for creating features. For example, darts can be 
introduced into a T- spline by inserting a few adjacent rows of control points with zero knot 
intervals, as shown in FIG. 15. Having two adjacent knot intervals of value zero introduces a 
local triple knot, and the surface becomes locally C° at that knot. The shape of the crease is 
5 controlled by the placement of the inserted control points. The crease can be made less sharp 
by replacing the zero knot intervals in FIG. 15 with small non-zero knot intervals. 

Converting a T-spline into a B-spline surface 

This refinement algorithm makes it easy to convert a T-spline in Si into an equivalent 
10 B-spline surface in S n by simply computing the transformation matrix Mi jn as discussed in the 
section on T-spline spaces. 

A standard T-spline is one for which, if all weights Wj=l , then 

n n 

J^sWjBifat) = J* M B£s 9 i) si. This means that the denominator in equation (2) is identically 

equal to one. Hence, the basis functions provide a partition of unity and the T-spline is 
15 polynomial. Thus, an algebraic statement of necessary and sufficient conditions for a T-spline 
to be standard is each row of Mi )n sums to 1 . 

The insertion algorithm can produce a surprising result: a T-spline for which not all 

n 

Wj=l but yet ^ w,5,O,0 = 1 . This T-spline can be called a semi-standard T-spline. FIG. 16 

shows two simple examples of semi-standard T-splines. The integers (1 and 2) next to some 
20 edges are knot intervals. A semi-standard T-spline space S is one for which there exists some 

n 

elements of S for which w i B i (s, f ) s 1 » and not all wj=l . A non-standard T-spline space is 

m 

n 

one for which no elements exist for which ^ w i B l (s,t) s 1 . These definitions are more 

precise because they allow for the notion of a rational (weights not all =1) T-spline that is 
either standard, semi-standard, or non-standard. The distinction is made based on which type 

25 of T-spline space it belongs to. 

FIG. 17 illustrates an overview of an embodiment of a method for defining bicubic 
spline surfaces that provides local refinement to control meshes in a computing environment. 
One operation that is included in the method is specifying knot intervals associated with the 
spline control mesh as in block 250. Another operation is imposing a local knot coordinate 

30 system based on the knot intervals as in block 252. The local knot vectors for the control 
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points may be inferred in order to produce basis functions for the control points as in block 
254. A further operation is inserting a single control point into the control mesh without 
altering the surface as in block 256. This insertion of the control point can permit partial 
rows of control points to be inserted that terminate in a T-junction which is a special control 
5 point. 

Extracting Bezier Patches 

It is advantageous to represent in Bezier form the patches that comprise a T-spline, 
10 because a tessellation method may then be applied with the minor modification that each T- 

junction must map to a triangle vertex in the tessellation to assure that cracks will not appear 

in the tessellation. In general, there are fewer Bezier patches in a T-spline than in the 

equivalent B-spline surface. 

The domains of the Bezier patches that comprise a standard T-spline can be 
15 determined by extending all T-junctions by two bays, as illustrated in FIG. 18. The rectangles 

in FIG. 18b are Bezier domains. The reason for this can be understood by considering the 

knot vectors for the basis functions of each control point. 

Bezier control points can be obtained by performing repeated local knot insertion. 

Recall that a B-spline surface can be expressed in Bezier form using multiple knots, and that 
20 a zero-knot interval implies a double knot. For the knot interval configuration in FIG. 19b, 

the 4X4 grid of control points surrounding F are the Bezier control points of that patch. Thus, 

the Bezier control points for face F in FIG. 19b can be determined by performing local knot 

insertion. 

25 

Merging B-splines into a T-spline 

This section discusses how to merge two B-spline surfaces with different knot vectors 
into a single T-spline. Often in geometric modeling, portions of an object are modeled 
independently with different B-spline surfaces that have different knot vectors, such as the 
30 hand in FIG. 2. FIG. 20 illustrates the problem. The control grids are defined over different 
knot vectors. Merging them into a single B-spline requires that they have the same common 
knot vector, so knot insertion must first be performed before merging can proceed. As FIG. 
20c illustrates, however, those required knot insertions can significantly increase the number 
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of rows of control points. If additional surfaces are subsequently merged onto these two 
surfaces, the number of knot lines can multiply further. 

One prior solution to the problem of merging two NURBS surfaces into a single 
surface without a proliferation of control points was to use cubic NURSSes introduced by 
5 Thomas W. Sederberg in 1998. Since cubic NURSSes allow different knot intervals on 

opposing edges of a face, two NURBS control grids can be merged into a single control grid 
without propagating knot lines. FIG. 21 shows the result of merging two identical B-spline 
cylinders with different knot vectors. Unfortunately, any NURSS representation introduces an 
unsightly bump at the junction of the two cylinders. This failed attempt at solving the merge 

10 problem using NURSSes helped motivate the creation of T-splines. 

The procedure using T-splines is illustrated in FIG. 22. For a C" merge (n € {-1, 0, 1, 
2}), n + 1 columns of control points on one patch will correspond to n + 1 columns of control 
points on the other patch. We consider first the C° merge in FIG. 22a. To begin with, each B- 
spline has triple knots (double knot intervals) along the shared boundary, as shown. For a C° 

15 merge, one column of control points will be shared after the merge. If the knot intervals for 
the two T-splines differ along that common column, control points are located along the 
boundary edge so that the knot intervals agree. In this example, the knot intervals on the red 
B-spline are 1, 3, 2, 2 and on the blue B-spline are 2, 1, 4, 1. After inserting offset control 
points on each control grid along the soon-to-be-joined columns as shown, the common 

20 column of control points has knot intervals 1,1,1,1,2,1,1. 

Typically in this process, the control points that are to be merged will have slightly 
different Cartesian coordinates. For example, A on the red patch might differ slightly from A 
on the blue patch. Simply take the average of those two positions to determine the position of 
the merged control point. 

25 A C 2 merge is illustrated in FIG. 22b. The basic idea is the same as for a C° merge. 

The differences are that four knot intervals a, b, c, d must correspond between the two 
surfaces, as shown. Also, three columns of control points must be merged, as shown. FIG. 23 
shows the results of a C° and a C 1 merge. FIG. 4 shows an application of this merge 
capability in a NURBS model of a human hand. This is a merge. 

30 



T-NURCCs 

T-NURCCs are NURCCs (Non-Uniform Rational Catmull-Clark Surfaces) with T- 
junctions in their control grids and using T-spline concepts. NURCCs are generalizations of 
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tensor product non-uniform B-spline surfaces: if there are no extraordinary points, if all faces 
are four-sided, and if the knot intervals on opposing edges of each face are equal, NURCCs 
degenerate to non-uniform B-spline surfaces. NURCCs are also generalizations of Catmull- 
Clark surfaces: if all knot intervals are equal, then NURCCs degenerate to a Catmull-Clark 
5 surface. Likewise, T-NURCCs are a generalization of T-splines, Catmull-Clark surfaces, and 
cubic NURBS surfaces. 

NURCCs enforce the constraint that opposing edges of each four-sided face have the 
same knot interval and NURCCs have stationary refinement rules. It is also for this reason 
that NURCCs are capable of local refinement. 

10 The refinement rules for NURCCs are thus identical to the refinement rules for 

NURSSes if we require opposing edges of each four-sided face tQ have the same knot 
interval. Those refinement rules are discussed in Thomas W. Sederberg in 1998. 

We now discuss how T-junctions can be used to perform local refinement in the 
neighborhood of an extraordinary point. To simplify our discussion, all extraordinary points 

1 5 are separated by at least four faces, and all faces are four-sided. These criteria can be met by 
performing a few global refinement steps, if needed. Thereafter, all refinement can be 
performed locally. For example, any suitably large regular sub-grid of a NURCC control grid 
can undergo local knot insertion, as discussed previously. Also, refinement near an 
extraordinary point can be confined to the neighborhood of the extraordinary point. 

20 To explain how to perform local refinement in the neighborhood of an extraordinary 

point, a way to perform local knot insertion in the neighborhood of a single (valence four) 
vertex in a T-spline will be explained. Referring to FIG. 24, the refinement begins with the 
black control grid. Then, it is legal using procedures defined previously to insert all of the 
control points in row 1, followed by row 2, followed by column 3, followed by column 4. 

25 Then the additional control points can legally be inserted in like order, etc. What is produced 
is a local refinement in the immediate neighborhood of one central control point. Note that 
this refinement scheme can split faces at any ratio p . For a valence-4 point, changing p does 
not change the limit surface since we are merely doing B-spline knot insertion, but when we 
adapt this scheme to extraordinary points, p will serve as a shape parameter. FIG. 25 shows 

30 the effects of changing p . 

We now present the local refinement rules for T-NURCCs at an isolated extraordinary 
point. Referring to FIG. 26, knot interval di is split into knot intervals d\ is split into knot 
intervals p d\%xv& (1 - p)d\\ likewise for the other knot intervals adjacent to the extraordinary 
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point. If p = 1/2 and if all the initial knot intervals are equal, the limit surface obtained using 
this local refinement is equivalent to a Catmull-Clark surface. 

Lower-case letters refer to knot intervals and upper-case letters to points. Vertices A, 
B, C, D, Q, R, S, T are the initial control points, prior to refinement. After refinement, these 
5 vertex points are replaced by new vertex points denoted with primes: A', B', C", D', Q', R' . 



[e 0 + (1 - P)e x ] [{d 0 + (1 - p)d x )A + (pd i +d 2 )B] 
(d 0 +d l +d 2 )(e 0 +e ] +e 2 ) 

| [pe x + e 2 ) [(d 0 + (1 - p)d x )C + (pdj +d 2 )D] 
(i/ 0 +^+f/ 2 )(e 0 +e,+e 2 ) 



1 0 There are three types of edge points: E, H, and G. 



E 2 =pM 2 +(\-p) 




where 
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2pd l +d 2 +h A 



B + 



2d 0 +2(\-p)d i 



2(d 0 +d l ) + d 2 +h 4 



2{d 0 +d 1 )+d 2 +h 4 




| K, +d Q +(l~ P)d x ] [(pg, + e 2 )D + (gg + (1 - p)e, )B] 
(</_, +d 0 +d x )(e 0 +e,+e 2 ) 
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Edge point G, = 



e 0 +e x +e 2 



pg. +e 2 



R + 



g 0 +(l-/g)g| 
e Q +e t +e 2 
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There are five different types of vertex points: those that replace A, B, Q, D, and R. 
We will denote the new vertex point at A by A' , etc. 



A' = p 2 A + 2p(l-p) 



1=0 



E 
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n-1 



where n is the valence, m, = (/im + /i,.2 + /in-2)/2, and _/J = hj-yh^- 
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5'=(l-p) 



+ P 



M „ , d -i +d o+Q-P) d \ B 



d_ x + d 0 +d x 



d_ x +d 0 +d x 



, [pd x R + (d_ x +d 0 +(l-p)d x )D][e_ x +e 0 +(l-p)e,] 



R'=- 



pe x 



e_ x +e 0 + e, 



(d_ x +d 0 +d x )(e_ x +e 0 +e x ) 
o , g-i+g 0 +0-/>)g| ^ 



Local refinement at such an extraordinary point is illustrated in FIG. 6, which shows a 
T-NURCC that has undergone four steps of local refinement. The yellow dots highlight four 
1 5 T-junctions. Note that this locally-refined mesh has two orders of magnitude fewer faces than 
it would have using global Catmull-Clark refinement. 

This discussion has assumed that extraordinary vertices are separated by at least four 
faces, and this can be accomplished by performing a few preliminary global refinement steps. 
It is possible to derive local refinement rules that would not require such initial global 
20 refinement steps, but there are additional special cases to consider. 

Away from extraordinary points, NURCCs are C 2 , except that zero knot intervals will 
lower the continuity. At extraordinary points where all edges have the same knot interval 
value, an eigenanalysis for the valence 3 case shows X x =\>^, 2 =X z >^ i >... where, for 
example, 
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l + 2p-h > 9 2 -hV(p 2 + 6p + l)(p^l) T 

This gives analytic proof that for the valence three case, the surface is G ] for any legal value 
of p ( 0 < p <1). A similar result can be obtained for valence five. For higher valence, the 
symbolic algebra expressions get unwieldy, but a sampling of specific values of p has not 
5 revealed any case where the G 1 condition is not satisfied. If the knot intervals for edges 
neighboring on extraordinary point are not equal, then the empirical evidence suggests G 1 . 

T-splines and T-NURCCs permit true local refinement: control points can be added 
without altering the shape of the surface, and (unless there are knot intervals with a value of 
zero) the new control points can be moved and the surface will remain C 2 . Since T-NURCCs 

10 generalize NURBS and Catmull-Clark surfaces, a modeling program based on T-NURCCs 
can handle any NURBS or Catmull-Clark model as a special case. 

The present invention has been illustrated using bi-cubic B-spline surfaces. However, 
the present invention can also be implemented using B-spline surfaces of any degree. The 
concepts explained in this invention would enable one skilled in the art to derive T-splines of 

15 any degree, and this is a straightforward exercise. One difference between even-degree T- 
splines and odd-degree T-splines is that for even degree T-splines, the knot intervals are 
associated with control points and for odd-degree, the knot intervals are associated with edges 
of the control grid. 

An additional use of the present invention is to use T-NURCCS to provide a seamless 
20 representation for the intersection of two T-NURCCs, or of any special case of T-NURCCS, 
such as bicubic NURBS surfaces or T-splines. 

It is to be understood that the above-referenced arrangements are illustrative of the 
application for the principles of the present invention. It will be apparent to those of ordinary 
skill in the art that numerous modifications can be made without departing from the 
25 principles and concepts of the invention as set forth in the claims. 



